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STABLE  EQUILIBRIA  IN  A  SCALAR  PARABOLIC  EQUATION 
WITH  VARIABLE  DIFFUSION 


by 

G.  Fusco  and  J.  K.  Hale 


Abstract 


A  scalar  parabolic  equation  with  nonconstant  diffusion  and  nonlinear 
source  term  is  considered  and  some  aspects  of  the  influence  of  chang 
ing  the  diffusion  on  existence,  stability  and  bifurcation  properties 
of  the  equilibria  are  discussed. 


1 .  Introduction 


We  deal  with  existence  stability  and  bifurcation  properties  of  equi¬ 
libria  of  the  problem 


(1) 


ut  “  (cVx  +  f(u)*  x  6  t-1’1) 

ux(-l,t)  =  ux(l,t)  =  0 


where  c  >  0  is  a  continuous  function  and  f  is  C*. 

The  initial  value  problem  for  (1)  is  well  posed  in  the  Sobolev  space 
H*(-l,l),  [1],  and  any  bounded  orbit  approaches  an  equilibrium  as  t  -*■  *>  ([3], 

[5], [7]).  Therefore  a  basic  problem  in  understanding  the  dynamic  of  (1)  is 
the  description  of  the  set  of  equilibria  of  (1)  and  of  the  way  this  set 
changes  with  the  diffusion  function  c  and  with  the  source  term  f.  Related 
important  problems  are  the  characterization  of  the  pairs  (c,f)  such  that 
(1)  has  stable  nonconstant  equilibria  and  to  understand  the  role  of  bifurca¬ 
tion  in  the  appearance  of  stable  equilibria. 

For  any  nonlinear  function  f,  Chafee  [6]  proved  that  when  c  is  con¬ 
stant,  no  stable  nonconstant  equilibrium  exists.  Chafee's  result  was  generaliz 

ed  by  Hale  and  Chipot  [2]  that  showed  that  the  same  result  holds  true  if 
2 

c  €  C  and  c^  <_  0.  Finally,  Yanagida  [10]  has  shown  that  if  c  is  written 
2 

as  c  *  a  ,  a  >  0  a  necessary  and  sufficient  condition  for  the  nonexistence 
of  a  function  f  such  that  (1)  has  a  stable  nonconstant  equilibrium  is  that 
axx  <,  0.  Other  results  concerning  the  existence  of  stable  nonconstant  equili¬ 
bria  are  due  to  Matano  ([7], [9])  that  has  shown  that,  if  f  is  a  cubic  poly- 

3 

nomial  as  f  *  u  -  u  and  c(x)  >_  1  on  intervals  [-!,<*], [6,1]  and  <_  e 


on  [y»$]  a  <  y  <  6  <  6  and  e  is  sufficiently  small,  then  (1)  has  a 
stable  nonconstant  equilibrium.  Fife  and  Peletier  [13]  have  also  considered 
equations  related  to  (1)  which  have  stable  nonconstant  equilibria. 

For  the  n  dimensional  version  of  problem  (1)  in  a  bounded  domain  n 
and  with  constant  diffusion,  Casten  and  Holland  [11]  and  Matano  [8]  have 
shown  that, if  Q  is  convex, any  stable  equilibrium  must  be  a  constant. 

3 

Matano  has  also  shown  that,  assuming  f  of  the  type  f  =  u  -  u  ,  for  some 
nonconvex  domains,  there  exist  stable  nonconstant  equilibria.  Hale  and  Vegas 
[4]  have  shown  the  existence  of  stable  nonconstant  equilibria  for  a  large 
class  of  nonlinearities  and  for  domains  £1£  that  can  be  considered  as  pertur¬ 
bations  of  a  domain  which  is  the  union  of  two  disjoint  convex  domains. 

We  assume  that  c  is  even  and  f  is  odd  and  such  that 


(2) 


f (0)  =  f(l)  =  0 
f (u)  >0  for  u  €  (0,1) 
f(u)  <  0  for  u  €  (l,00), 
f'(0)  ^  0,  f'(l)  *  0. 


(see  fig.  1.) 

Under  these  assumptions  we  give  an  estimate  of  the  number  of  equilibria 
of  (1)  in  terms  of  c  and  f .  We  prove  that  for  any  f  of  type  (2) ,  if  c 
is  sufficiently  close  to  the  step  function 


(3) 


c 


1  for  x  €  [-1,1]  U  [1,1], 

Cq  >  0  for  x  €  (-1,1),  0  <  l  <  1 


and  Cq  is  sufficiently  small,  problem  (1)  has  at  least  a  pair  of  stable  non¬ 
constant  odd  monotone  equilibria.  Finally,  we  show  that  if  c  =  c^  depends  on 


a  parameter  y  €  [0,1]  and  is  an  equilibrium  of  (1)  with  exactly  k 
zeros  that  bifurcates  at  y  =  0  from  the  zero  equilibrium  and  becomes  stable 
at  y  *  1  then  as  y  goes  from  0  to  1,  u^  must  go  through  at  least  k 
secondary  bifurcations. 


2.  Existence 


We  are  interested  in  studying  problem  (1)  for  c  in  the  set  of 
continuous  and  positive  functions  c:  [-1,1]  F.  Nevertheless,  for  the 
analysis  that  follows,  in  particular  for  the  discussion  of  stability  where 
we  consider  function  c  that  are  "close"  to  the  step  function  c  defined 
by  (3),  it  is  convenient  to  study  problem  (1)  for  a  wider  class  of  diffusion 
functions  c  that  are  allowed  to  have  jump  discontinuities  at  a  finite 
number  of  fixed  points  in  [-1,1].  To  keep  the  notation  simple  and  since  we 
suppose  c  even,  we  consider  only  the  case  of  two  points  of  discontinuity 
at  x  =  ±4,  0  <  4  <  1.  Everything  we  say  extends  to  the  case  of  any  number 
of  points  of  discontinuity. 

Let  if  be  the  set  of  nonnegative  even  functions  c:  [-1,1]  -*■  R  which 
have  continuous  restrictions  to  [0,4)  and  to  [4,1]  and  possess  the  left 
limit  c(4  )  of  c(x)  as  x  -*■  4.  For  any  c  €  if  let  J^cF  be  the  set 

J£  d=f  {(x,y)|x  =  4,  y  €  [c(4~),c(4)]}  and  C  d=f  Jc  U  graph  c.  We  sup¬ 
pose  that  if  is  endowed  with  the  topology  associated  with  the  following 
notion  of  convergence  that  allows  a  sequence  of  continuous  functions  to  con¬ 
verge  to  a  function  that  has  a  jump  at  x  =  ±4:  we  say  that  cn  €  if, 
n  =  1,...,  converges  to  c  G  if  if  and  only  if  Cn  converges  to  C  in  the 

sense  of  Hausdorff  as  n  <*>. 

The  class  of  diffusion  function  that  we  are  going  to  consider  is  the 
subset  if+  cif  defined  by  the  condition  inf  c  >  0.  Clearly,  if  is  a  dense 
subset  of  X  and  if  we  assume  in  ^  the  topology  of  uniform  convergence  in 

[-1,1],  then  5f  is  continuously  embedded  in  if+.  Henceforth  we  allow  c  in 


problem  (1)  to  be  a  generic  c  €5s  .  This  requires  that  (1)  be  complemented 
with  the  jump  conditions 


c(±*+)ux(±t+)  =  c(±£')ux(±0, 

therefore  the  equilibrium  problem  corresponding  to  (1)  becomes 


(4) 


(cux)x  +  f(u)  =  0 
ux(-i}  *  ux(l)  =  0 
c(±4+)ux(±£+)  =  c(±£")ux(±£~). 


and  reduces  to  the  standard  problem  for  c  €  5^\ 

By  letting  u  =  u,  v  *  cux>  problem  (4)  transforms  in  the  equivalent 

system 


(5) 


u  = 

X  c' 


vx  -  -f(u), 


vC-1)  •  v(l)  •  0. 


Note  that  the  jump  conditions  express  just  continuity  of  v  at  x  =  ±£ 
and  therefore  they  are  included  in  the  requirement  that  u,v  be  continuous 
in  [-1,1]. 

The  hypothesis  on  f  and  a  maximum  principle  argument  imply  that 
solutions  of  (4)  or  (5)  satisfy  -1  <_  u(x)  <_  1,  therefore  we  can  also 
assume  that  f  is  bounded  so  that  the  solution  u(c,a,x),  v(c,a,x)  of  the 
initial  value  problem 


I  u(-lj  =  a,  v(-l)  =  0 


is  defined  for  all  (c,a,x)  e  #  *  [-1,1]  *  [-1,1], 


Lemma  1.  u(c,a,x),  v(c,a,x)  are  continuous  functions  of  (c,a)  uniforml 


m  x  and  possess  a  continuous  first  derivative  with  respect  to  a  and 


also  with  respect  to  x  except  possibly  at  x  =  : 


The  proof  of  this  lemma  is  a  standard  application  of  the  general  theory 
of  differential  equations. 


To  discuss  the  existence  of  space  dependent  equilibria  of  (1),  i.e., 
the  existence  of  nonconstant  solutions  of  (4),  we  note  that  these  solutions 
are  in  one  to  one  correspondence  with  the  a  t  -1,0,1  such  that  v(c,a,l)  * 
If,  for  a  i  0,  we  let  6(c,a,x)  be  the  angle  (positive  clockwise  around 
x  in  space  u,v,x)  swept  by  the  vector  u(c,a,x')  defined  by 


”u(c,a,x') 

u(c,a,x’)  =  v(c,a,x')  , 


when  x*  goes  from  -1  to  x,  then  a  necessary  and  sufficient  condition 
in  order  that  v(c,a,l)  *  0  for  some  a  f  -1,0,1  is  that  5(c,a,l)  be 
equal  to  irk  for  some  integer  k  i  0. 

The  angle  6(c,a,x)  can  be  defined  also  for  a  =  0  so  that  <s(c,a,x) 
is  continuous  in  (c,a,x).  In  fact  by  performing  the  polar  coordinate  trans 


formation  u  =  p  cos  6,  v  =  -p  sin  6  it  is  found  that  6(c,a,')  is  the 


solution  of  the  problem 


(7) 


c(x) 


•  2r 
sin  <5 


+  f(u(c,a,x)) 
u(c,a,x) 


2. 

cos  0, 


fi(-l)  =  0. 


Moreover,  since  by  lemma  1,  u(c,a,x)  is  continuous  in  (c,a)  uniformly 
in  x  and  u(c,o,x)  =  0,  it  follows  that  if  (c',a)  converges  to  (c,o) 
in  *  [-1,1],  f (u(c' ,a,x))/u(c' ,a,x)  converges  uniformly  to  f ' (0)  in 

[-1,1].  This  implies  that  as  (c',a)  -*•  (c,0),  6(c',a,x)  converges 

uniformly  to  the  solution  6(c,0,*)  of  the  problem 


(8) 


sin26  +  f'(0)cos26, 

X  c(x) 

<5(-l)  =  0. 


From  (7), (8)  and  lemma  1  it  also  follows  that  6(c,a,x)  is  continuously 
differentiable  with  respect  to  a.  We  also  note  that  6(c,±l,x)  =  0  and 
that,  for  a  €  (1,1),  6(c,a,x)  is  an  increasing  function  of  x  because 
the  right  hand  sides  of  (7) , (8)  are  >  0. 

For  later  use  we  also  introduce  the  angle  o(c,a,x)  which  is  defined 
as  6(a,x)  with  the  vector  u(c,a,x)  replaced  by  its  derivative  u  (c,a,x) 
with  respect  to  a.  It  may  be  useful  to  note  that,  if  E  is  the  surface 
in  the  space  of  u,v,x  defined  by  the  solutions  of  (6),  then  u  (c,a,x)  is 
tangent  to  the  cross  section  of  E  at  x  at  the  point  (u(c,a,x),  v(c,a,x),x) 
(see  fig.  2). 


It  is  easily  seen  that  ofCja,-)  is  the  solution  of  the  problem 


(9) 


12  2 
°x  =  F(x7  sin  0  +  (u(c,a,x))cos  a, 

a(-l)  =  0. 


Now  consider  an  interval  [-t,S.]c[-l,l]  and  let  c  be  the  supremum 
of  c(x)  in  [-&,£].  Then  for  x  €  [-£,£]  (8)  implies 

(10)  6x  >  — -  sin26  +  f'(0)cos26. 

c 

Since  a  simple  computation  shows  that  the  solutions  of  (10)  with  the  equality 

- 1/2 

sign  increase  of  tt  each  time  that  x  increases  of  tt  (f  *  (0) / c)  ,  from 

(10)  and  the  fact  that  6(c,a,x)  is  a  nondecreasing  function  of  x,  it 
follows  that 

6(c,a,l)  >_  tt.  Integer  part  of 

This  estimate  together  with  the  continuity  of  6  (c  ,  • ,  1 )  and  the  fact  that 


5(c,±l,l)  =  0  imply 


Theorem  1.  The  number  N  of  nonconstant  equilibria  of  (1)  satisfies  the 
condition 


(ID 


N  >  2. 


Integer  part  of 


Remark .  In  the  proof  of  theorem  1.  no  use  was  made  of  the  eveness  of  c 

and  oddness  of  f.  Thus  theorem  1.  holds  for  generic  c,f.  We  also  note 

that  the  conclusion  of  theorem  1.  is  also  true  if  [-£,£]  is  replaced  by  a 

measurable  set  Ec  [-1,1]  of  measure  2£. 

Let  sfc  =  (a|6(c,a,l)  =  kir).  The  set  s^  can  be  identified  with  the 

set  of  equilibria  of  (1)  that  have  exactly  k  zeros.  If  the  right  hand  side 

of  (12)  is  >_  2k,  then  s^  is  nonempty  and  by  mean  of  equation  (7)  it  is 

possible  to  obtain  some  information  on  the  "shape"  of  equilibria.  To  this  aim 

let  0  <  u  <1  and  a  €  s^  be  given  and  f  c  [-1,1]  be  the  set  where 

|u(a,x)|  <  u.  To  get  a  bound  for  the  measure  of  J  n  (-£,£)  we  let  (xi»x^) 

be  the  smallest  interval  containing  ^  fl  (-£  £),  L  =  x2  -  Xj  its  length,  and 
—  def 

n  =  min_f(u)/u,  then  by  applying  to  (7)  the  same  procedure  used  for  deriving 
|u|<u 

(10)  from  (8)  we  obtain 


i  2  _  2  _ 

6  >_  —sin  6  +  n  cos  6  ,  x  €_/ D  (-£,£) 

x  c 


and  therefore,  since  6(c,a,x)  is  an  increasing  function  of  x. 


4 


S(c,a,x2) 


d6 


2—2 
6(c,a,x.)  _l_sin  6+ncos  6 


-10- 


From  this  estimate  it  follows 


(12) 


L  <  kir 


n 

because  S(c,a,X2)  -  S(c,a,xp  <  kir.  The  estimate  (12)  shows  that 
inf{ |u(c,a,x)  | }  converges  in  measure  to  1  in  (-£,£)  as  c  +  0. 


For  k  /  1  nothing  can  be  said  on  the  behavior  of  solutions  in  s^  outside 
the  interval  (-£,£)  and  solutions  in  s^  could  be  almost  trivial  in  the 
sense  that  they  could  be  near  zero  outside  (-T,T)  and  oscillate  in  (-£,£). 
This  can  not  happen  when  k  =  1  because  solutions  in  s^  are  monotone  and 
therefore  if  there  is  a  point  x  €  (-T,£)  where  |u(c,a,x)|  is  near  1  the 
same  is  true  in  [-l,x]  or  in  [x, 1]. 

In  what  follows  we  are  interested  in  solutions  of  (4)  that  are  odd 
functions  of  x.  It  is  easily  seen  that,  due  to  the  assumption  that  c  is 
even  and  f  is  odd,  when  on  the  basis  of  (11)  it  is  possible  to  conclude 
that  Sj  is  nonempty,  then  it  also  contains  at  least  a  pair  of  odd  solutions 
that  transform  in  each  other  under  the  transformation  x  -x.  Clearly,  if 
u(c,a,*)  is  one  of  these  odd  solutions,  and  |u(c,a,x)+l|  <  e  in  [-l,x] 
then  |u(c,a,x)±l|  <  e  in  [-jc,l],  Therefore  on  the  basis  of  (12)  we  have 


Theorem  2.  For  any  c  such  that  the  right  hand  side  of  (11)  is  %  2,  problem 
(1)  has  an  equilibrium  which  is  an  odd  and  increasing  function  of  x.  If  c 
is  deformed  so  that  c  0,  then  all  odd  increasing  equilibria  of  (1)  converge 
to  the  function 

[  -1  for  x  €  [-1,0) 

0  for 


z 


1 


for 


x  =  0 
x  €  (0,1] 


-11- 


and  the  convergence  is  uniform  in  compact  sets  in  [-1,0)  U  (0,1]. 

In  the  statement  of  theorem  2.  and  in  the  following,  we  always  refer 
to  the  increasing  equilibrium,  with  it  being  understood  that  there  is  also 
a  decreasing  equilibrium  that  transforms  into  the  other  one  under  the  trans 
formation  x  -x. 

3.  Stability. 

Let  c°e£  be  the  function  defined  by 

1  x  €  [-1,-1)  U  (1,1], 

~0 
c  = 

0  x  €  (-!,£) . 

, 

In  this  section  we  prove  the  following 

Theorem  3.  Let  f  be  a  continuously  differentiable  odd  function  that 
satisfies  (2) .  Then  there  is  a  set  W  ci  if  such  that 
(i)  W  is  open  and  connected  in  if 
(ii)  c°  belongs  to  the  closure  of  W  in  if 
(iii)  for  any  c  6  W  problem  (1)  has  an  odd  increasing  (and  an  odd 
decreasing)  equilibrium  which  is  stable. 

Note  that  theorem  3.  implies 

Corollary.  For  any  odd  C^-function  f  that  satisfies  (2)  there  is  a 
c  €  if  such  that  problem  (1)  has  a  stable  nonconstant  equilibrium. 


To  prove  theorem  3.  we  need  a  few  lemmas. 


Proof.  If  S/  is  empty,  the  lemma  is  obvious.  Therefore  we  assume  that 
is  nonempty.  Then  there  exist  c  €  ig',  a  €  (-1,1)  {0}  and  k  >  0  such 

that 

6(c,a,l)  =  kn, 

and  the  largest  eigenvalue  \  of  problem  (13)  is  negative.  If  one  lets 
w  *  r  cos  v,  cw^  =  -r  sin  v  in  equation  (13),  it  is  found  that  r, v  satisfy 

(14)  rx  =  sin  vcos  v  (f '  (u(c,a,x))  -  -  X)  r, 

(15)  vx  =  sin2v+  (f*  (u(c,a,x))-X)cos2v, 


with  the  boundary  conditions  v(-l)  =  0,  v(l)  =  in  for  some  integer  i. 


But  i  must  be  zero  because  the  eigenfunction  w  corresponding  to  the 
largest  eigenvalue  never  vanishes  and  therefore  v  must  stay  in  the  interval 
(-ir/2,ir/2) .  Since  X  is  negative,  for  v  =  a,  the  right  hand  side  of  (9)  is 
always  smaller  than  the  right  hand  side  of  (15);  therefore,  we  have 
cr(c,a,l)  <  v(l)  =  0.  On  the  other  hand,  equation  (9)  implies  o(c,a,x)  >  - 
It  follows  that  sin(6(c,a,l)  -  o(c,a,l))  i  0.  Since  the  derivative  of  6 
with  respect  to  a  is  related  to  a  by 

P  <Sa  2 

(16)  — 2~ 2 — J  =  Sin  (6"a) 

P  Sa+Pa 

it  results  6  (c,a,l)  i  0.  Then  the  lemma  follows  by  the  implicit  function 

cl 

theorem. 

Remark .  Lemma  3  is  actually  a  special  case  of  a  general  situation.  In  fact, 
the  largest  eigenvalue  being  negative  for  an  equilibrium  point  uQ,  implies 
the  semigroup  generated  by  the  linear  variational  equation  is  exponentially 
asymptotically  stable.  Thus,  small  perturbation  in  c  will  yield  another 
equilibrium  point  near  Uq  also  has  the  largest  eigenvalue  negative  and  it 
will  be  stable. 

In  the  proof  of  lemma  3,  we  have  seen  that  a(c,a,l)  <0  is  a  necessary 
condition  for  X  to  be  negative.  We  note  that  this  condition  is  also  suf¬ 
ficient.  This  follows  from  the  fact  the  solution  of  (15)  depends  continuously 
on  X,  coincides  with  o(c,a,-)  for  X  =  0*  and  increases  unboundedly  as 
X  -►  for  x  j*  -1.  Therefore,  if  a(c,a,l)  is  negative,  there  exists  a 
unique  negative  X^  such  that  the  solution  of  (15)  vanishes  at  1.  Then, 
if  r(*)  is  any  non  zero  solution  of  (14),  w( • )  =  r(-)cosv(*)  is  an  eigen- 


K»|=) 
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function  of  (13)  that  does  not  vanish  in  [-1,1].  Thus  Xq  <  0  is  the 
largest  eigenvalue  of  (13).  Therefore,  we  can  state 


Proposition  1.  A  necessary  and  sufficient  condition  in  order  that  the 


largest  eigenvalue  of  problem  (13)  be  negative  is  that  o(c,a,l)  be 
negative. 

Lenma  4.  Let  c  €.*  be  a  function  of  type  (3)  and  (1).  be  problem  (1) 


with  c  »  c.  Then  if  c„  is  sufficiently  small,  prob 


1).  has  an  o 


corresponding  linear  problem  (13)-  is  negative. 

4i2 

Proof.  By  theorem  1,  if  cQ  <  f '  (0)  there  exists  an  a  €  (-1,1)  n,  (0) 

such  that  u(c,a,»)  is  an  increasing  equilibrium  of  (1)-.  The  same  condition 

c 

together  with  the  evenness  of  c  and  the  oddness  of  f  ensure  that  a  can 

also  be  chosen  so  that  u(c,a, • )  is  an  odd  function.  To  prove  that  the 

largest  eigenvalue  of  the  linearized  problem  at  u(c,a, •)  is  negative  if  cQ 

is  sufficiently  small,  we  recall  [12]  that  the  eigenvalues  of  (13)  do  not 

decrease  if  f' (u(c,a,x))  is  replaced  by  a  function  q(x)  >_  f' (u(c,a,x)) .  It 

follows  that,  if  we  let  q  *  max_  f' (u(c,a,£)) ,  it  suffices  to  show  that 

u>u(c,a,A) 

for  Cq  small,  the  largest  eigenvalue  of 


Wxx  +  <*w  *  Xw  »  x  e  C-1*-*-)  U  (t,l) 
c0wxx  +  f’ (u(c,a,x))w  *  Xw,  x  6  (-t,t) 

wx(-l)  »  wx(l)  *  0 

cowx^'A+^  "  V-O. 

\(t*>  ■  Vx(,‘> 


is  negative. 


,**  p‘«  ,%  *. 


From  a  result  of  Yanagida  [10],  it  follows  that  the  largest  eigen¬ 
value  of  this  problem  is  negative  if  there  is  a  strictly  positive  function 


Wq  that  makes  the  left  hand  sides  of  (17)  equal  to  zero,  satisfies  the  last 
two  equations  (18)  and  moreover,  is  such  that 

(19)  woxC_1)  <  °»  w0x(1^  >  °* 


We  look  for  an  even  such  w_  and  therefore,  we  assume  w  (0)  =  0  and 

0  ox 

consider  only  the  interval  [0,1].  Since  f ' (1)  <  0  and  by  theorem  2, 
u(c,a,£)  ■+  1  as  cQ  -*■  0,  q  is  negative  for  small  value  of  cQ.  Therefore, 
if  wQ  exists,  in  the  interval.  [i,l],  it  must  have  the  expression 

(20)  wQ(x)  =  A  sinh  [(-q)1/2(x-£)]  +  B  cosh  [(-q)1/2(x-t)]. 


and  the  coefficients  A,B  must  satisfy  the  conditions 

(21)  A  >  0,  |  >  -  tanh  [(-q)1/2(l-£)] 


ensuring  that  wQ(x)  is  positive  in  [£,1]  and  wox(l)  >  0. 

To  compute  wQ(x)  in  the  interval  [0,£],  we  must  solve  the  problem 


(22) 


c0Woxx  +  f* (u(c,a,x))wQ  =0,  x  €  (0,£), 
wQ(0)  =  C,  w0x(0)  =  0, 


where  C  is  a  positive  constant  to  be  chosen  later. 


From  now  on  we  set  for  simplicity  u  =  u(c,a,*)»  u  =  u(£),  u  =  u(l).  In  order 
to  solve  this  problem  we  must  overcome  the  difficulty  lying  in  the  fact  that 
u  is  only  known  to  be  an  odd  increasing  solution  of  problem  (4)-.  To  this 


end  we  observe  that  since  u  is  increasing,  we  can  perform  the  change  of 

variable  x  =  u  *(u)  5(u).  By  making  this  change  of  variable  in  (4)~ 

c 

and  by  observing  that  the  oddness  of  u  implies  5(0)  *  0,  we  see  that 
5  satisfies 


-co  ^  +  f<u> 


u  G  (0,u) , 


+  f(u)  = 


u  6  (0,u), 


5(0)  =  0 


5(u)  =  i. 


lim  5(u)  »  1,  lim  £'(u)  =  00 . 
u*5  ih-u 


By  using  the  fact  that  ^  one  sees  that  the  same  change  of  variables 


applied  to  (22)  yields 


-p "o  ■  cop "o  *  f'“  ■  “•  u  6  (0-u>- 
«0(0)  -  c,  ujco)  >  0, 


where  wQ  has  been  identified  with  the  function  wQ(5(*)).  From  (23)  it 
follows  that,  for  u  €  (0,u) 


C  1 1 

=  2  f  (s)ds  +  K  d=f  g (u) , 

5  '  >\i 

f1 

where  K  >  -2  f(s)ds  is  an  integration  constant  that  together  with  u,u 


satisfies  the  conditions 


f(s)ds) 


*  11 


cQg(u)  =  2l_  f  (s)ds. 


which  correspond  to  the  last  two  conditions  (24).  Sine  g*  =  -2f  and  the 

t 

first  equation  (23)  implies  that  the  coefficient  of  w^  in  equation  (25) 
is  equal  to  -f,  equation  (25)  becomes 


ii  • 


gw0  ‘  fw0  +  f'W0  =  (gW0  +  fw0}  =  °* 


Thus  gwQ  ♦  fwQ  =  const  =  0  because  wQ(0)  =  0  and  f(0)  *  0. 

It  follows  that,  with  a  proper  choice  of  the  constant  C  appearing 
in  (25), 


1/2 

w0  =  8  ’ 


(for  u€  [0,u]). 


If  the  expressions  (20) (27)  are  patched  together  at  x  =  l  (corresponding 
to  u  «  u)  by  imposing  the  conditions 


it  is  found  that 


w0(*+)  -  w0(O.  wQx(£+)  =  c0w0x(£'). 


A  =  (g(u))1/2  ,  B  =  -c01/2  £(^TT2  ■ 

(-q)  ' 


Therefore,  it  follows  that  if 

c01/2f(u)  _  1/2 

(28)  ■ "  -  ,  H - — T7T  <  tanh  [(-q)4/^(l-£)] , 


then  wQ  satisfies  all  the  conditions  ensuring  that  the  largest  eigenvalue 
of  problem  (17) (18)  is  negative.  We  shall  prove  that  this  is  the  case  for 
Cq  sufficiently  small.  The  proof  is  a  discussion  of  the  asymptotic  depen¬ 
dence  of  u,g(u)  on  cQ  defined  by  equations  (26)  for  cQ  ->■  0. 

By  the  change  of  variables  u  =  u  +  (u-u)t,  s  =  u  +  (u-u)o,  the 
first  of  equations  (26)  transforms  as 


(29) 


2l/Z  JO  VO  fi(x,o) 


(x-a+  (l-x)o)dcr 


\  —  1/ 2 
J  dx 


1  -  t. 


where 

6(x,o)  =  -[(1-u)  +  (u-u) (1-x) (1-a)] , 

a  =  >  I. 

u-u 


As  u,u,  also  6  and  a  depend  on  c„.  Let  a  =  lim  sup  a.  Since 

u  c  0 

_  o  _ 

u  -*>  1,  u  ■+  1  as  cQ  -*•  0,  the  above  expression  of  6(x,a)  implies  that 
6(x,o)  -*■  0  uniformly  as  c^  -►  0.  Therefore,  f (1)  =  0  implies  that  the 
ratio  f (l+6(x,a))/6(x,a)  converges  uniformly  to  f* (1)  /  0.  This  and 
equation  (29)  imply  o  <  «. 

We  also  have 


lim  inf  - 3-j  [_  f(s)ds  =  lim  inf  [_  = 

Cq  -►  0  (1-u)  J  u  cQ  -►  0  a  J\x  (u-u) 


(30) 


with  7*  (t)  «  (u-1)  ♦  (u-u)t.  Since  7'  (x)  -*>  o  uniformly  as  cQ  -*•  0, 


the  ratio  f(l+T' (t))/?' (t)  -*>  f*(i)  <  o  uniformly  as  cQ  -*■  0.  From  this 


and  the  fact  that  a  is  >  1,  it  follows 


(31) 


lim  inf  ^  f  (T-a)dT  >  -f'(l)f  (l-x)dT  «  -  ±f'(l) 

crt  -*■  0  J 0  6'(t)  J0 


Equations  (30) (31)  and  the  second  condition  (26)  imply 

crtg(u) 


(32) 


lim  inf 


cQ  -  0  (l-u)‘ 


-ZT2  1  -f’d) 


Therefore,  by  taking  into  account  that  lim  q  «  f'(l),  we  obtain 


C<T° 


co1/2f(u) 


Urn  sup  ■  _  jy2 - 3172  *  - *  i/2  lim  SUP  c0  -  i/S'l  i/2 

cn-0  (-q)1/2(g(u))1/2  (-f'(l))1/Z  cQ+0  0  (1-u)  c01/2(g(u))1/2 


that  is:  the  left  hand  side  of  (28)  converges  to  zero  as  Cq  0.  This 
proves  the  lemma  because  the  limit  of  the  right  hand  side  is  >0. 

Proof  of  theorem  3.  By  lemma  4,  there  is  a  number  e  >  0  such  that,  if 

Y  c  is  the  curve  y  (b|0  <  cQ  <  e)  and  c  €  y,  then  problem  (1)~ 

has  an  odd  increasing  equilibrium  which  is  stable.  Since  by  lemma  3  Sf  is 

open  in  5f+  there  exists  an  open  neighborhood  W  of  y  in  5f+  such  that 

for  c  €  W  problem  (1)  has  a  stable  equilibrium  u  .  It  is  easy  to  see  that 

c 

W  can  be  chosen  so  that  uc  is  odd  and  increasing.  In  fact  from  the  proof 
of  lemma  3  it  follows  that  for  c  in  a  neighborhood  of  c  €  y,  u£  is  the 
only  equilibrium  in  a  neighborhood  of  u» .  On  the  other  hand  the  evenness 
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of  c  and  the  oddness  of  f  imply  that  also  v  defined  by  v  (x)  = 

c  c 

-u  (-x)  is  an  equilibrium  of  (1).  Since  u  is  odd  and  u  -*•  u  as 
c  c  c  c 

c  -*•  c  also  v  converges  to  u„  as  c  -»■  c.  This  contradicts  uniqueness 
c  c 

of  u  unless  u  is  odd  and  therefore  proves  oddness.  Since  u  is 
c  c  c 

close  to  u  it  vanishes  only  at  x  =  0.  From  this  and  the  fact  that 
c 

solutions  of  (4)  with  only  one  zero  are  monotone,  it  follows  that  u£  is 
increasing.  The  mapping  Cq  -*•  c  is  continuous  as  a  map  from  (0,e)  into 
Sf+.  Therefore  y  is  locally  compact  as  a  subset  of  if+ .  Thus,  by  standard 
arguments  there  is  a  continuous  function  $  :  (0,1)  -►  jf  such  that  the  curve 
y  -  (c|’c  =  <|»(s),  s  €  (0,1)}  is  contained  in  W  and  c^  is  in  the  closure 
of  y  in  if ’.  Since  W  is  open  and  if  is  continuously  embedded  in  , 
w  nif  is  open  in  Sf.  From  this  and  the  continuity  of  <J>  it  follows  that 
there  is  a  subset  W  c  W  fl  Sf  which  is  open  and  connected  in  if  and  contains 
y.  Since  c°  is  in  the  closure  of  y  in  jf  the  proof  is  completed. 


4.  Secondary  bifurcation 

In  this  section  we  consider  a  family  (c  )  c  if  of  diffusion  functions 


ue[-i,i] 


def 


c^  depending  continuously  on  a  parameter  u.  We  let  6(y,a,x)  =  6(c^,a,x) 

and  assume  that  c  <  c  for  y_  >  y  and  that  6 (0,0,1)  =  kir  for  some 
U2  U1  1  L 

k  >  0.  Then  equation  (8)  implies  that  6(y,0,l)  <  kit  for  y  <  0  and 
6(y,0,l)  >  kir.  If  we  also  assume  that  f  satisfies  the  condition 


(33) 


f (u)  <  f»(0)u 


U  £  (0,1], 


then  from  (7)(8)  it  follows  6(y,a,l)  <  6(y,0,l)  for  a  +  0,  y  €  [-1,1]. 
Therefore  s^  is  empty  for  y  <_  0,  nonempty  for  y  >  0.  Thus  y  =  0  is  a 
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bifurcation  point.  It  is  easy  to  see  that  in  this  situation,  for  y >  0 
and  small,  s^  contains  solutions  that  are  small  and  converge  to  zero  as 
y+0,  i.e.,  solutions  that  bifurcate  from  the  zero  solution.  These  solu¬ 
tions  are  unstable  for  y  small  because  the  largest  eigenvalue  of  problem 
(13)with  u  e  0  is  f ' (0)  >  0  and  the  eigenvalues  of  (13)  are  continuous 
functions  of  c£^  On  the  other  hand  we  have  seen  in  theorem  3  that  if 
c  is  suitably  chosen  then  there  exist  stable  nonconstant  equilibria  of 
(1).  Therefore,  it  can  be  .expected  that,  if  u^  is  a  continuous  function 
of  y  €  [0,1]  such  that  u^  =  0,  u^  is  a  solution  of  (4)  in  s^  for 
y  €  (0,1],  and  Uj  is  stable,  some  kind  of  secondary  bifurcation  takes 
place  at  somey  €  (0,1).  This  conjecture  is  true,  we  have  in  fact  the 
following. 

Theorem  4.  Suppose  that  cy  is  as  before,  f  satisfies  (33) ,  ug  is  an 

equilibrium  of  (1)  which  is  equal  to  zero  for  y  =  0,  has  exactly  k 
cy 

zeros  for  y€  (0,1],  depends  continuously  on  y  and  u^  is  stable  (in  the 
sense  that  the  largest  eigenvalue  of  the  linearized  problem  at  Uj  is 
negative) ,  then  there  exist  numbers  0  <  y^  <...<  y^  <  1  such  that  y^, 
i  =  l,...,k  is  a  bifurcation  point. 

def 

Proof.  Let  a^  =  u^(-l).  Then  =  0  and  therefore  (7), (8)  imply 

a(0,ag,l)  =  6(0,ag,l).  We  also  have  6(y,ay,l)  =  kiT  for  y  €  (0,1],  thus 
by  the  continuity  of  6  with  respect  to  c,a  and  the  continuity  of  c^, 
uy  with  respect  to  y,  it  follows  o(0,ag,l)  =  kir.  On  the  other  hand 
Proposition  1  and  the  stability  of  Uj  imply  o(l,aj,l)  <  0.  Therefore, 
by  continuity  there  exist  0  <  y^  <  . . .<  y^  <  1  such  that 


'22- 


o(y  ,a  ,1)  =  (k-i)it  i  =  1, . . .  ,k. 

1  yi 

Moreover,  it  is  obvious  that  y^,  i  *  can  be  chosen  so  that  in 

any  neighborhood  of  y.  there  ekist  y”  <  y.  <  y  such  that  o(y,a  ,1)  >_ 

1  1  y 

(k-i)ir  >  o(y, a_,l).  This  on  the  basis  of  the  geometrical  meaning  of  the 
y 

angle  a  implies  that  y^  is  a  bifurcation  point. 

Theorem  4  says  that  going  through  k  secondary  bifurcations  is  a 
necessary  condition  in  order  that  an  equilibrium  with  k  zeros  that  bi¬ 
furcates  from  the  zero  solution  becomes  stable.  From  the  proof  of  the 
theorem  and  Proposition  1  it  follows  that  if,  as  y  goes  from  0  to  1, 
u^  experiences  exactly  k  bifurcations  at  0  <  yj  <  . . .<  y^  <  1  each  one 
of  which  is  simple  in  the  sense  that  at  any  yj^  two  new  solutions  bifurcat¬ 
ing  from  appear,  then  Uj  is  stable  (see  fig.  3  for  che  case  k  =  2) . 

This  observation  shows  that  in  a  certain  sense  the  converse  of  theorem  4  is 


also  true. 
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